Hypersingular integral equations with applications 
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NUMERICAL TREATMENT OF HYPERSINGULAR INTEGRAL 
EQUATIONS WITH APPLICATION TO CRACK MECHANICS 
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In this paper we present a treatment of hypersingular integrals and integral equations. It 
is shown the equivalence of different definitions of hypersingular integrals. Then we construct 
explicit analytical solution to a characteristic equation and prove convergence of the colloca- 
tion method to this solution. Further, we expand these results to full hypersingular kernels, 
and give a demonstration of the developed method to some problems of crack mechanics 
including crack problem for porous elastic materials, in the framework of Cowin-Nunziato 
model. 

1. The method. Let us consider the hypersingular equation of the following 
type 
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{x~ty 



Ko{x,t) 



dt = f{x) 



\x\ < 1 



where Ko{x,t) is a regular part of the kernel. Direct numerical treatment of Eq.(l.l) 

is not easy. A particular case is the hypersingular equation with the characteristic 
kernel 

' 9{t)dt 



xe{a,b), f (x) e C2 (a, b) 



;i.2) 



We find a solution of this equation bounded on the both ends x = a,b. The latter can 
be defined in exact explicit form, which is given by the following inversion formula 



[x — a) {b — x) 
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f{t)dt 
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This result shows that any bounded solution of Eq.(1.2) vanishes at x ^ a,b. To 
construct a direct collocation technique to numerically solve equation (1.2) for arbitrary 
right-hand side, we divide the interval (a, b) to n small equal subintervals by the nodes 
a = to, ti, t2, ...,tn-itn = b. Then the length of each interval is h = {b — a) /n, and 
the nodes tj = a + jh, j = 0, l...,n. Let us denote the central points of each sub- 
interval {ti-i,ti) by Xi, so that Xi = a + {i — 1/2) h, ^ = 1, ...,n. Therefore, we try to 
approximate Eq.(1.2) by the linear algebraic system 
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It is proved in [1] that by assuming x G (a, b), the difference between solution g (x) of 
the system (1.4) at the point x and solution given by explicit formula (1.3) tends to 
zero, when n —>■ oo. It can be proved that system (1.4) admits the explicit solution 
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If the number of nodes increases then g{xi) — > g{x). 
Let us we consider the full equation 



{x - ty 



+ Ko {x, t) 



where 



Ko(x,t) = 



g (t) dt ^f(x), xe (a, b) 
dKi (x, t) 



dx 
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(1.6) 
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Its bounded solution can be constructed by applying inversion of the characteristic 
part, that reduces Eq.(1.6) to a second-kind Predholm integral equation: 

b 

g{x) + j N, {x, t) g (t) dt = fi{x), xe (a, b) , (1.8) 



where 



Ni{x,t) = 
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I ^(r-a) (b-T) (x-t) 



Ki (r, t) dr 
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y — a) {b — t) {x — t) 
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It is known from the classical theory of the Cauchy-type integrals that if / (x) G 
Ci (a, b) , Ki {x, t) e Ci [(a, b) x (a, b)] then also /i (x) e Ci (a, b) ; A^i (x, t) e 
Ci [{a, b) X (a, b)]. Further, we prove that if / {x) G Ci (a, 6); Ki {x, t) G Ci [(a, b) x (a, b)], 
then for any x G (a, b) the difference between solution g {x) of the linear algebraic sys- 
tem 
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Xi tj Xi tj — 1 
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9 {tj) = f (xi) , i 
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and the bounded solution of equation (1.6) tends to zero when /i — > (i.e. n — > oo). 

2. Application to crack mechanics. The linear Cowin-Nunziato theory of 
homogeneous and isotropic elastic material with voids is described by the following 
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system of partial differential equations (where (f) — u — uq is the change in volume 
fraction from the reference one [2,3]) 

fjL A u + {X + n) grad div u + P grad = 

(2.1) 

aA(j) — ^(p — f3 div u = , 

where fi and A are classical elastic constants; a, j3 and - some constants related to 
porosity of the medium. Besides, u denotes the displacement vector. The components 
of the stress tensor are defined, in terms of the functions u and 0, by the following 
relations {5ij is the Kronecker's delta) 

aij = A 5ij SfcA: + 2 /i Eij + (5 (f) Sij 

If we consider a plane-strain boundary value problem for the thin crack of the length 
2a with plane faces, dislocated over the segment —a<x<a along the x-axis. Let the 
plane-strain deformation of this crack be caused by a normal load — ctq symmetrically 
applied to the faces of the crack. For the last problem the boundary conditions over 
the line y — are 



xy^O, ^ = {\x\ < oo) , ayy = -(7o {\x\ < a) , Uy ^ {\x\ > a) , (2.3) 



dy 
where 

with the first two numbers c , H being dimensionless and the quantities Zi , I2 - of 
dimension of length. 

Let us apply the Fourier transform along the x-axis to relations (2.1)-(2.3). Then 
the problem can be reduced to the following integral equation: 



j gmix - m =-{^-Nf^, \x\<h. (2.5) 
-h 

where g is the opening of the crack face, N = il^/li) H , {0 < N < 1) and 
1 °° I I 

K{x) = — [ -p-[2Nc's\q-\s\) + {l-N){l-N-c^)q]e~'''-'ds 
ZTT J q{s) 
—00 

1 °° 

= - / L{s) cos{sx)ds, q = q{s) = Vs^ + l-N , ^^'^^ 

TT J 


L(s) = [2Nch\q - s) + (1 - N){1 -N- c^)q 

The kernel (2.6) admits explicit representation by special functions that permits direct 
estimate of its singular properties. It can be seen that the kernel is hypersingular, and 
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we apply the proposed method to numerically solve equation (2.5). Then we study in 
detail the influence of the porosity of the material to the stress concentration coefficient. 
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